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In this paper, we analyze the stability of a homogeneous self-gravitating plasma, having a non-zero resistivity. 
This study provides a generalization of the Jeans paradigm for determining the critical scale above which gravitational 
C • ' collapse is allowed. 

We start by discussing the stability of an ideal self-gravitating plasma embedded in a constant magnetic field. We 
outline the existence of an anisotropic feature of the gravitational collapse. In fact, while in the plane orthogonal 
["t"! . to the magnetic field the Jeans length is enhanced by the contribution of the magnetic pressure, outside this plane 
perturbations are governed by the usual Jeans criterium. The anisotropic collapse of a density contrast is sketched in 
details, suggesting that the linear evolution provides anisotropic initial conditions for the non-linear stage, where this 
effect could be strongly enforced. 

The same problem is then faced in the presence of non-zero resistivity and the conditions for the gravitational 
collapse are correspondingly extended. The relevant feature emerging in this resistive scenario is the cancellation of 
(— > . the collapse anisotropy in weakly conducting plasmas. In this case, the instability of a self-gravitating resistive plasma 
is characterized by the standard isotropic Jeans length in any directions. The limit of very small resistivity coefficient 
is finally addressed, elucidating how reminiscence of the collapse anisotropy can be found in the different value of the 
perturbation frequency inside and outside the plane orthogonal to the magnetic field. 
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1. Introduction 

In many astrophysical and cosmological systems the presence of a plasma component has a very important role in 
determining the shape and the behavior of their equilibrium configurations. The peculiarity of these plasma configura- 
tions with respect to those ones observed in laboratory, relies on the dominant character of the gravitational interaction 
in determining the stability properties. Indeed, as firstly suggested by Jeans [1, 2], the gravitational interaction is able 
to induce the collapse, as long as a critical scale of the configuration is reached. Such a scale, depending on the sound 
speed and on the mass density of the medium is commonly known as the Jeans length. On the other hand, dealing 
with a self-gravitating plasma instead of a fluid brings in the equilibrium all the typical features observed in magneti- 
cally confined and highly ionized gases, like the emergence of Alfven and magnetosonic waves. Moreover, we stress 
that the possibility to postulate the presence of a magnetic field, is allowed since it is observationally demonstrated 
by the direct observation of astrophysical systems (see e.g. Refs. [3, 4, 5, 6, 7, 8, 9, 10] and, more recently, Refs. 
[11, 12, 13, 14, 15, 16, 17, 18]; see also Refs. [19, 20, 21, 22, 23] for review works). 
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The fact that, in our analysis, both the background mass density and magnetic field of the configuration are taken 
homogeneous is justified by the often slow variation of these quantities in real astrophysical systems, even over scales 
for which the self-gravity is already relevant (for instance, the primordial cosmological plasma and the ionized inter- 
galactic baryonic component [24]). Indeed, our study concerns the linear stability of a homogeneous magnetized and 
self-gravitating plasma, endowed with a finite value of the resistivity coefficient. This latter dissipative feature is here 
introduced to account for the non-ideal nature of the most commonly observed space plasma. Significant reconnection 
processes of the magnetic profile are often observed or argued via the interpretation of data from astrophysical con- 
figurations. Despite the effects of a finite resistivity coefficients are particularly important in the non-linear regimes, 
where the establishment of a turbulent profile of the plasma can phenomenologically enforces the resistivity (see for 
instance the question concerning the so-called anomalous resistivity in the configurations of stellar accretion disks 
[25] from the plasma instabilities raised from the streaming of electrons, and also [26, 27, 28]), nonetheless, we will 
show how its presence is crucial already in the linear case, when dealing with the stability properties. 

As a first step, we analyze the linear stability for the ideal case, when the resistivity coefficient of the plasma 
vanishes. In this limit, we essentially reconstruct the Jeans paradigm for the gravitational stability of the plasma 
structure. The stability out of the plane orthogonal to the constant magnetic field remains still characterized by the 
same Jeans length obtained originally for the fluid scheme. A relevant new feature emerges, however, in the plane 
perpendicular to the magnetic force straight lines, where the contribution due to the magnetic pressure affects the 
equilibrium enhancing the value of the Jeans length by a term corresponding to the square of the Alfven velocity 
in the plasma. Such additional contribution enters the Jeans length expression on the same footing as the sound 
speed contribution and therefore its relevance strictly depends on the ratio of the sound speed to the Alfven one. 
In particular, the greater the Alfven speed is, the larger is the anisotropy in the gravitational collapse, inside and 
outside the orthogonal plane. We properly describe this effect by following the behavior of a over-dense region 
during the linear evolution, which is accordingly squeezed on the orthogonal plane. Indeed, in the linear regime, the 
density contrasts grows without a real gravitational collapse (that takes place essentially in the non-linear stage of 
the evolution), but the growth is slower on the plane where the magnetic pressure affects the Jeans scale. It is worth 
noting that, in the perturbation scheme, the magnetic pressure is, despite its name, anisotropic, being provided by 
the scalar product between the background and the perturbed magnetic field. It is just this feature that introduces the 
anisotropy in the perturbation evolution. The mode that becomes unstable on the orthogonal plane corresponds to the 
quasi-stationary one (typical of the slow magnetosonic configuration of a non-gravitating plasma). The presence of 
gravity alters the nature of this mode, making the system unstable, but with a greater Jeans length with respect to the 
directions out of this plane, for which the quasi-stationary mode would be absent in the non-gravitational case too. 

The gravitational stability analysis is then faced taking into account a non-vanishing resistivity coefficient. In 
this case, our work outlines how the stability condition is now fixed by the request that the scale length should be 
smaller than the Jeans value, independently from the direction along which the perturbation mode propagates, as in 
the standard model. This issue relies on the plasma dispersion relation in presence of a non-zero resistivity. However, 
the anisotropic evolution of the linear perturbations is not completely removed, since their growth proceeds at a 
different velocity that depends on the angle between the perturbation wave number and the background magnetic 
field. This feature is well elucidated in the limit of very small (but non zero) values of the resistivity coefficient, i.e., 
for a very conducting fluid. In fact, although the orthogonal stability is not present when an arbitrarily small value of 
the resistivity is taken onto account, the time evolution of the perturbations is still affected by a certain dependence on 
the direction, as far as the resistivity coefficient is not too large. This is due to the small, though real, values that the 
frequency takes on the orthogonal plane to the magnetic field when the resistivity is very small, in comparison to the 
other directions. 

It is important to remark that in the resistive MHD model proposed in this paper, any other dissipative effects 
have been neglected (like those driven by viscosity) although in many astrophysical contexts these effects are clearly 
relevant. Indeed, a correct and more realistic picture of plasma instability mechanism should be properly given by 
the visco-resistive MHD approach. In fact, the relative magnitude of the viscous and magnetic diffusion rates can be 
parameterized through the magnetic Prandtl number. For the cosmological plasmas, for instance, the viscous diffusion 
is more important than resistive diffusion, at recombination and for more recent times. Thus, viscosity dominates 
on resistivity at all scales of cosmological interest and becomes relevant (i.e., the associated frequency becomes 
equal to the Alfven frequency) at recombination for scales of the order of 10 Kpc comoving or less. Contrarily, 
the corresponding scale for the resistivity is very small and it is under the scales relevant for cosmology. However 
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both these "critical" scales depend on time and, therefore, the viscosity always dominates over resistivity, but both 
dissipative effects can indeed be neglected when addressing the propagation of Alfven waves [29]. 

Summarizing, our study outlines the role that the plasma nature of the system plays in the gravitational instability 
dynamics. We generalize the concept of Jeans length and outline the anisotropic dynamics induced by the magnetic 
pressure, both in the ideal and in the resistive case. The obtained results acquire particular interest when they are 
treated as initial conditions for the non-linear and turbulent collapsing dynamics. The paper is organized as follows. 
In Sec. 2, we provide the fundamental equations for the resistive MHD equilibrium and develop the perturbation 
scheme. The resulting differential system is combined in order to obtain a single ordinary differential equation for the 
fluid velocity, which gives the dispersion relation describing the system stability. Then, we discuss the zero-resistivity 
case in Sec. 3 and the case with finite resistivity Sec. 4. In both cases we analyze the evolution of the density contrast in 
momentum and coordinate space, in order to outline the growth of the anisotropy. Finally, discussion and conclusions 
follow in Sec. 6. 

2. Resistive perturbative MHD dynamics 

2.1. Basic equations 

In this Section we recall the basic equations of resistive MHD. The mass conservation and the configuration of the 
Newtonian gravitational field are described by the continuity and Poisson equations 

d ( p + V-pv = 0, V 2 - 4-nGp = , (1) 

respectively, where p is the mass density, v is the velocity field, <t> is the gravitational potential and G is Newton 
constant. Moreover, the single-fluid dynamics of the plasma is described by the Euler equation in the presence of a 
magnetic field B, i.e., 

pd,v + p(v ■ V)v + VP + pVd> - (V X B) X B/4tt = . (2) 
The electromagnetic interaction is summarized by the Maxwell equations for the electric (E) and magnetic fields 

d t B + c(VxE) = 0, (3) 
4^-c(VxB) = , (4) 
V ■ B = , (5) 

and by the Ohm law for the density current J 

B + *-**=V. (6) 

c 

Combining Eqs. (3), (4) and (6) we get the fundamental relation for the evolution of the magnetic field: 

d,B = V x (v x B) + fjV 2 B , (7) 

where we have introduced the diffusion coefficient f\ = rjc 2 /4n. For our analysis, the dynamics is summarized by 
Eqs.(l), (2) and (7). 

2.2. Perturbation scheme 

In order to analyze the implications that the physics of a resistive plasma has on the structure formation mechanism, 
we now follow the standard perturbation approach, as in the Jeans model [30]. In this respect, we consider a small 
perturbation around a static and uniform solution (indicated by the subscript 0). Thus, we assume p = po + p\, with 
pi <K po (the same is also valid for P, v, <1> and B) where the background is characterized by constant po, Po, Q>o, 
Bo and vo = [24, 31]. We follow the standard procedure that takes this as a solution of the above equations even 
if the Poisson equation is not actually satisfied (unless po = 0). This fact, that corresponds to neglect the effects of 
gravitation in the unperturbed solution, is known as the "Jeans swindle" (see for example [24]). 
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In this picture, the system describing the first-order perturbative dynamics reads now 



d,pi +p V-Vi =0, (8) 

V 2 *! - \nGp x = , (9) 

podtVi + v 2 V Pl +poV*i - (V x Bi) x B Q /(4n) = , (10) 

0,1*1 - f r l 2 B l + B (V ■ vi) - (Bo ■ V)vi = , (11) 

VBi=0, (12) 

where the pressure and density perturbations are related through the adiabatic sound speed, i.e., Pi = v 2 pi . Introducing 
the dimensionless magnetic fluctuation b = B\/Bq (where Bo = Bo Bo), we consider perturbations iff of the following 
form 

ilKr,t) = ^(t)e- i,l T , (13) 

where \f/\ stands for [p\, V\, b, <I>i] and q — qq denotes the wave-number (in the following we will drop the tilde for 
the sake of simplicity). This way, the corresponding evolution equations in the harmonic space can be simply obtained 
by the substitution V — > iq. A further simplification can be obtained by decomposing V\ as Vi = v\q + v\ (where 
v\ ■ q — 0), and by introducing the following variables 

6 = Pl /p , 6 = i{q-vi), b = b-B , v = iq(v\-B Q ). (14) 

In this scheme, the perturbative dynamics results to be described by the following system which involves only 
scalar quantities 

(5 + (9 = 0, (15a) 
8 - u>l6 - aj 2 A b = , (15b) 
b + 6(1 - ft 2 ) - fiv + u~b = , (15c) 
i + IMolb = , (15d) 

where we have defined 

B 2 

fj. = B ■ q e [0, 1] , tol = v]q 2 -AnGpa , u A = -r—q = v A q , (o„ = fjq 2 . (16) 

47rp 

The system (15) can be now reduced to a unique equation for the variable 6 

6» (4) + fjq 2 6 {3) + (to 2 + u>l) <9 (2) + T)q 2 oj 2 <9 {1) + p 2 ^ 2 = 0, (17) 

where denotes the t th derivative with respect to time. Since (17) is a linear, fourth-order, ordinary differential 
equation with constant coefficients, the explicit solution is of the form 9(t) oc e la " . The same consideration stands for 
all the variables of Eq.(14) also. This way, we get the following algebraic equation for the variable <±> 

<jj 4 + ia^u? + a2(i) 2 + ia\u + «o = , (18) 



where we have set 

fl = P 2 ^1<^o > fl i = w >; w o ' fl 2 = - (w 2 + to 2 ) , «3 = -to n . (19) 

Since all the quantities admit the same time dependence, e.g., 6 oc e mt , the gravitational stability is determined by 
the sign of to 2 , i.e., if the imaginary part of to results to be negative. A solution with to 2 > (< 0) corresponds to 
an oscillating (exponential) perturbation, while for to — one gets a static density contrast. In the following, we will 
focus on the analysis of the solutions of Eq.(18) discussing in some details the different cases. 
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3. Zero-resistivity case: analytical and numerical analysis 



Let us start by considering the ideal MHD limit. When rj = 0, Eq.(18) admits the following solutions: 



a> 2 ± (q) = \[{u>l + (o 2 A ) ± A /( W 2 + W 2 ) 2_ 4/Z 2 W 2 W 2 j _ (20) 

We recall that oj^ can be positive or negative, depending on the relative strength of the pressure and gravitational force. 
In fact, the standard Jeans picture is recovered when a> 2 = 0: in this case, the solution is simply to 2 — a> 2 and the 
oscillating (or > 0) and collapsing (oj 2 < 0, taking the negative solution) modes are separated by the critical (Jeans) 
wave-number 



1J 



4nGp (21) 



In the presence of a magnetic field (a> 2 + 0), it can be shown that 

[ sign(w? + a?) fx = , 
signal) = J (22) 
1 sign(w ") p + . 

Here, the relation for p = is evaluated by considering the a>+ solution since, in this case, a> 2 results always positive 
yielding an oscillating density contrast evolution (and vice versa for^u + 0). This implies that, when the perturbation 
wave-vector lies in the plane orthogonal to the background magnetic field (p = 0), the magnetic pressure can induce 
stable perturbations that would otherwise collapse in the absence of a B-field. In fact, even if to 2 < (indicating 
an instability according to the standard Jeans criterion), nevertheless a> 2 + a> 2 > if co A is sufficiently large. The 
corresponding critical wave-number q c is readily found to be 



4nGpQ 



+ V? 



Ou = 0) . (23) 



In particular, perturbation modes with q > q c are stable, while those with q < q c grow exponentially. Since q c is 
always smaller than q j, the presence of the magnetic field increases the stability of the system. By other words, the 
critical wave-number is simply obtained by replacing, in the expression for the Jeans wave-number, the sound speed 
with the effective speed v^.. — v 2 + v 2 . On the other hand, as soon as the perturbation is not exactly perpendicular to 
the background field (// + 0), the stability of the perturbation is dictated by the standard Jeans criterion, i.e., q c = qj. 

Concluding this section, let us now introduce for completeness the standard Jeans length and the magnetic critical 
length defined as 

*J = A , *JM = \i t, , (24) 

respectively. We recall that the second quantity defines the stability in the plane orthogonal to the background magnetic 
field, i.e., p — 0. 



3.1. Time evolution of the density contrast (Fourier space) 

In what follows, we numerically integrate the system (15), for a> n = 0, in order to characterize the behavior of the 
density contrast. The initial conditions are set as 6(to) = So, v(to) = 0, b{to) = and S\ to = 0. In Figure 1, we show the 
time evolution for perturbations with different wave-number, in the directions orthogonal and parallel to the magnetic 
field, respectively. It is evident that fluctuations with Aj < A < A jm, which would result unstable according to the 
standard Jeans criterion, are stabilized in the direction orthogonal to the background magnetic field. 
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Figure 1: Plot of the ratio S(t)/So as function of the time f, for/* = (left panel) and fj 2 = 1 (right panel). 



3.2. Time evolution of the density contrast (coordinate space) 

From the previous analysis, one would expect the gravitational collapse to be anisotropic. In this sense, in order 
to outline the nature of the anisotropy induced by the magnetic field in the physical space, we now consider an initial 
Gaussian over-dense region in real space, i.e., 

S( r , t = 0) = 6 c e-'~ ll,j2 , (25) 

where the background magnetic field Bq is directed along the y-axis. The corresponding initial over-density in Fourier 
space is still Gaussian (in this Section, we restore the tilde in order to denote quantities in harmonic space): 

* = 0) = 7t4t2 f 6 ^ 1 = 0y' 9 "V* 3 r = ~5 C e-* l2&2 , (26) 
(Inyi 1 J 

and, evolving such over-density in time, i.e., 6(q, f), we are able to (numerically) transform back to real space getting 
the profile of the density contrast at a given time. Taking v 2 A = 10v; (A jm ~ 3A j) and <x = 2Aj (so that Aj < cr < A jm), 
in Figure 2, we show the contour line of the density contrast where 5 = S c /e - Q315 c at different times. As a result, 
we find that the density contrast isolines tend to "squeeze" (for increasing time) along the direction orthogonal to the 
background magnetic field. 

This issue shows how the presence of the constant magnetic field confining the plasma can play an important role 
for the perturbation evolution since the collapse in the real space favors the formation of squeezed structures. Indeed, 
our analysis is limited to the linear evolution and we cannot make prediction about the final fate of such a collapsing 
profile, but it is remarkable that the initial conditions for the non-linear phase of this process results to contain an 
intrinsic anisotropic feature. 



4. Resistive case 



Let us now discuss in some details the stability properties of the the pure resistive MHD case. By setting 77 + in 
Eq.(18), the four solutions for a> write 



X/O- 



Figure 2: Plot of the isolines where S = 6 c /e, in the plane [x/cr, y/o~] (z = const) at different times (light gray=earlier, dark gray=later). 

where we have defined the following constants: 



In what follows, we consider oscillating and growing (or decreasing) perturbations separately and, for the sake of 
simplicity, we split co into its real and imaginary parts x and y, 



where x and y are both real. 

At a first instance, it is possible to show that real (non-zero) solutions (27) exist only in the case q > qj {col 0) 
and for fi = 1 . Recalling that 6 ~ e' CM , real frequencies or > correspond to the oscillating modes and, in this case, 
the solution results to be to = ±a>o and does not depend on the coefficient uj n . From this result, one can argue that pure 
oscillating perturbations are admitted only in the plane parallel to the background magnetic field, i.e., fi = 1, and for 
q > qj. In this regime, the modes are not affected by the presence of resistivity. 

Let us now turn our attention on the case of growing or decreasing perturbations. For scales smaller than the Jeans 
length and for fi + 1, the imaginary part y of at is always different from zero and, as we will show in the next Section, 
only solutions with y > are admitted (the regime y < does not occur). This fact results in a damping of the density 
contrast for q > qj. On the other hand, in the case of larger scales, i.e., q < qj, there are no oscillating perturbations 
and we find at least one growing exponential solution with y < and x — 0. In this sense, we can conclude that, 
also in the presence of resistivity, the basic Jeans instability criterion is confirmed and the critical wave-number which 
discriminates the gravitational collapse results, in turn, to be q = qj. For a better comprehension of these results, we 
now discuss in some details the different specific cases. Since the condition to have a growth of the density contrast is 
y < 0, in the following we do not consider purely real solutions. 




(28) 



(29) 



u> = x + ly , 



(30) 



4.1. Below the Jeans length: q> qj 

Imaginary frequencies: u = iy. We consider at first the purely imaginary solution where u - iy. We are interested in 
particular in defining the sign and the expression of y. It is convenient to introduce here the following notation 



± 1 
>W S T ± 4V^ 3 



(w 2 + S 2 w 2 ) + 40> 



V u r,-f(w5+S 2 w2) + 4o 



2 V2 



and 



2 V2 



(31) 



(32) 



where the Lundquist number S = to A /ai n = vjr\q = \nvjr\qc characterizes the properties of the plasma. In fact, at 
fixed scale, high values of S indicate strongly conducting plasmas. Let us also set 



pi/3 2 1/3 r 
+ 



P = I.+ V(£) 2 - 4(H 3 , 



3(2)1/3 3 pi/3 ' 

9 9/9 9 9 9\ / 9 9 9\ 9 9 / 9 9\ / 9 9 9\"~ 

S = 21u> u>~ [u> + /d"S u>~\ - Uo + S~co n ) 9u>~to n ( 1 + 8/i~S J - 2 (w + 5'~w^J 
T = 3w 2 w 2 (V5 2 - l) + (w 2 + s 2 w 2 ) 2 . 



Finally we name 



yi =y. 



yi = y + , yi = y!, 



J4 = )£. 



(33) 

(34) 
(35) 

(36) 



In the regime q > qj, purely imaginary solutions (when these exist) of Eq.(18) belong to the set {ryi, iyi, iyj, iy^}. It 
can be verified that all the y's are always greater than zero. This implies that, for q > qj, purely imaginary solutions 
always describe perturbations that are exponentially damped away. 

Complex frequencies: u> = x + iy. In this section we focus on the most general case x + and y + 0. It turns out 
again that y is always positive, while x = ±x s where 



u 2 a u n - 2 (w 2 + S 2 a> 2 )y + 3u n y 2 - 4y 3 
0)„ -Ay 



A very interesting case is given by the value y = where x 2 = x 2 , i.e., 
1 



xl = 



16 



8w 2 + (85 2 - 3) oj 2 n ± 2 ^16(4 + 32w 2 (l - 2/u 2 )S 2 co 2 + (16S 4 - 16S 2 + 3)aJ* 
^„ 2 =^o = ^(45 2 -l)/4. 



(37) 



(38) 



(39) 



As soon as we are considering the case aj 2 > 0, from Eq. (39) we must restrict our analysis to the case 5 > 1/2. We 
can identify the following different regions for the parameters: 

- 3 V2/8 < S < V5/4 and < yr 2 < y 2 : it results x 2 = x\ ; 

- S > V5/4 and <fi 2 < yU 2 : it results x 2 - x\; 

- S > V5/4 and yu 2 < yu 2 < fi 2 : it results x 2 = x 2 ± ; 

- ft 2 — pr : it results x 2 — x\ ; 



where we have defined 



2 49,cj 2 + {3 -16S 2 ) col 2 u 4 n (3 - 16S 2 ) + 16 (co 2 +S 2 co 2 n f 

^ ^ 256co 2 S 2 ' ^'^ 64co 2 co 2 S 2 ' (40) 

We can conclude that, in this parameter domain, the only significant information concerns the behavior of the 
frequency characterizing the oscillations of the perturbed quantities and, therefore, we can rule it out from the study 
of the gravitational stability. 

4.2. Above the Jeans length: q < qj 

In this Section, we consider the case q < qj. We find that there are no purely oscillating perturbations (y = 0), i.e., 
that the time dependence of the perturbation always contains an exponential part. In particular, there is always one 
and only one exponentially growing mode with y — y\ < and x — and at least one exponentially damped mode 
with y > and x = 0. This scenario is explored in detail in the following. 

Imaginary frequencies: co = iy. For q < qj, there is always one (and only one) solution of Eq.(18) with imaginary 
part y — y\ < 0. There is also at least another purely imaginary solution of Eq.(18), having y > and belonging to 
the set {iy2, iyi, iy*}. The existence of other purely imaginary solutions will depend in general on the value of the 
parameters fi, co 2 and S . However, this result shows that, above the Jeans length, an unstable mode is always present. 

Complex frequencies: co = x + iy. We focus now on the case x + and y + 0. The frequency co — x-v iy admits the 
solutions x — +x s (for y + co n j\ > 0, decreasing) and y where 



J6(2)2/3y + 12(2)1/ V/3u _ 3 r 8(5 2 w 2 + w 2) _ 3^ y 2/3 

h s T + nTn ■ < 41 > 



with 



v = 6 + Vtf 2 - 4v 3 , (42) 
6 = 2co 6 + 6co 4 (3 + (I - l2//)S 2 )col + 3co 2 S 2 (-3+2S 2 + / ?(9-24S 2 ))co 4 n + 2S 6 co 6 TI , (43) 
v = S A co 4 +2S 2 co 2 n u 2 {6p 2 + \) + col[ul-3co 2 ^. (44) 

The solution V2 is founded to be always positive. Thus we can conclude that for q < qj, complex solutions of Eq. (18), 
when they exist, always correspond to oscillating decreasing modes. Together with the result for the purely imaginary 
case, this implies that above the Jeans length, there is always one and only one unstable mode, independently on the 
values of the others parameters. 



5. Limiting cases of small resistivity and anisotropy 

The comparison between the resistive and non-resistive case underlines a difference that may seem odd at first; 
namely, the presence of an even infinitesimally small resistivity seems to completely erase the stabilizing effect (in 
the orthogonal plane) of the magnetic field, since in the resistive case we recover the standard, non-magnetic Jeans 
criterion for gravitational instability. One could then conclude that the support effect of the magnetic field is never ac- 
tually realized in practice, since an infinitely conducting fluid is obviously just a mathematical idealization. However, 
an analysis of the limiting case of "very strong" magnetic field and "small" resistivity shows that this is not the case, 
since the timescale for collapse on the orthogonal plane is much larger than outside that same plane, where it recovers 
the standard value, only dictated by the action of gravity. The orthogonal timescale becomes actually infinite when 77 
is exactly 0, reproducing the result of the non-resistive case described in the previous section. 

We then consider the four solutions of Eq. (27) in the limit of strong magnetic field and small resistivity, i.e., 
co 2 » co 2 , and co 2 » \co 2 \. By other words, we take the Alfven timescale to be much shorter than the resistive diffusion 
and standard gravitational collapse timescales. Thus we write co 2 = -e 2 co 2 with e <K 1 (we choose the minus sign 
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because we are interested in the Jeans-unstable solution), and we perform a series expansion of the frequency with 
respect to the small parameters e and e' = w^/w fl = 1/5 . We focus on the only Jeans-unstable mode, whose frequency 
in this limit reads (we only keep leading terms): 



oj = 



2 



—ie e da — i- 



-iea>A — —1 



0" = O); 
0" = 1), 



(45) 



in the direction orthogonal and parallel to the magnetic field, respectively. The mode is clearly unstable (Im(w) < 0) 
in both directions. However, the timescale over which the over-density grows is very different in the two cases. In 
the parallel direction, the leading term is of first order in the expansion parameters (and we recover the standard Jeans 
result or = a) 2 ); in the orthogonal direction, it is instead of third order in the expansion parameters. In particular, the 
frequency in the latter case is suppressed by a factor e/S with respect to the former, and actually becomes zero (i.e., 
the mode is stable) for an infinitely conducting plasma (5 = oo). If the timescale over which orthogonal perturbations 
grow is much larger than the time interval over which the system is observed, these perturbations can be considered 
stable for all practical purposes. In any case, they grow much slower than their parallely-oriented counterparts, and 
this open the possibility to observe an anisotropic collapse like in the case of non-resistive plasmas. 

5.7. Time evolution of the density contrast 

Let us now analyze the time evolution of the density contrast exploring the anisotropy of the gravitational collapse 
for different values of the resistivity r\. In this sense, we consider an initial Gaussian over-density in real space, i.e., 
of the form (25), and a background magnetic field Bo directed along the y-axis (as assumed for the zero resistivity 
case). The corresponding initial over-density in Fourier space has still a Gaussian profile. We evolve the Fourier 
over-density in time, as shown in Section. (3.2), in order to get 6(q; t) and then (numerically) transform back to real 
space. Moreover, we set v 2 A = \Qv 2 (Ajm ~ 3/1/) and cr — 2Aj (so that Aj < cr < Ajm)- 

We consider three different values for the resistivity, corresponding to a Lundquist number at the Jeans scale 
defined by 5 \ q=qj = {10, 0.3, 0.07). In Figure 3, we show the equal density contours where 6 = Q.996 c . As a result, we 
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Figure 3: Equal density contours corresponding to S = 0.995 c at different times (light gray=earlier, dark gray=later). We take v A = lOv^ 
(Ajm ~ 3/ij) and cr = 2Aj. From left to right, the Lundquist number at the Jeans scale is S(q = qj) = {30, 0.03, 1.5 X 10 3 [. 



show how the bigger the resistivity is the lower is the "squeeze" of the equal density contours along the orthogonal 
direction. In this respect, we conclude noting how the resistivity reduces the anisotropy in the density perturbations 
due to the presence of the magnetic field. 



6. Discussion and Conclusions 

We have considered the stability of a static, homogeneous infinite plasma. The perturbation equations and the 
stability problem has been studied for this system in some details and the perturbation frequencies have been classified, 
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pointing out different regions which are marked by the values of the resistivity and Alfven velocity. Two distinct 
regimes have been addressed: the ideal case and a resistive MHD scheme. 

In the ideal MHD picture, it has been outline how, if the propagation vector lies on the plane orthogonal to the 
background magnetic field, the stability is affected by the presence of the magnetic pressure, and, since we obtain 
v 2 s — > v 2 s + v 2 a , the collapsing characteristic scale is greater than the standard Jeans one. On the other hand, out of this 
plane, the standard Jeans criterion remains valid. The anisotropy in the density-perturbation profile has been shown 
as a consequence of the presence of the magnetic field: as time goes by, the equal density contours squeeze along the 
direction orthogonal to the background magnetic field. 

For the non-vanishing resistivity regime, the situation can be summarized as follows. The standard Jeans criterion 
governs the stability of all perturbations, irregardless of their direction. In fact, for small resistivity, the magnetic field 
lines are frozen into the fluid and the field-perpendicular motion will therefore drag the lines with it, and increase 
the magnetic pressure, while the field-parallel compression does not increase the magnetic field strength. However, 
in the limit of very small but non-zero resistivity, it is found that the speed of the perturbation growth is strongly 
direction-dependent. In particular, perturbations in the orthogonal plane grow much slowly with respect to those 
parallel to the magnetic field. This again leads to an anisotropy in the gravitational collapse similar to that observed in 
the zero-resistivity case. On the other hand, the difference in the evolution of the density perturbation profile along the 
orthogonal direction, as shown in the zero and small resistivity cases, disappears as the resistivity increases, suggesting 
a balance effect between the magnetic field and the dissipative term driven by resistivity. 

In this work, the resistive MHD system has been addressed as a simplified (one-fluid viscosity-free MHD) model 
describing the astrophysical plasmas. Indeed, one should complete this framework providing a generalization towards 
a more realistic models, e.g., visco-resistive MHD. Such scheme, allows to include the dissipation processes driven by 
viscosity and constitutes the basis of accretion disk models, since the (turbulent) visco-resistive plasma configurations 
can be regard as one of principal element for the momentum transport across the disk [25]. Moreover, a two-fluids 
description (with ions and electrons) or multi-fluids models with charged and neutral particles should account for 
other effects like the ambipolar diffusion, important for many astrophysical processes that involve lightly ionized gas. 
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